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Abstract 

Current-induced torque is formulated based on the spin continuity equation. The formulation 
does not rely on the assumption of separation of local spin and charge degrees of freedom, in contrast 
to approaches based on the s-d model or mean-field approximation of itinerant ferromagnetism. 
This new method would be thus useful for the estimation of torques in actual materials by first- 
principles calculations. As an example, the formalism is applied to the adiabatic limit of the s-d 
model in order to obtain the analytical expression for torques and corresponding (3 terms arising 
from spin relaxation due to spin-flip scattering and spin-orbit interaction. 
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I. INTRODUCTION 



Spin transfer torque is a torque acting on local spins as a result of an applied current. 
Such a torque has been discussed mostly based on s-d type of exchange interactioni^iii^Ji 
after the pioneering works by Berger-^ and Slonczewski^. In s-d models, the conduction 
electrons and localized spins are discriminated, and therefore the transfer of spin angular 
momentum between those two degrees of freedom occurs. However, in reality, this separa- 
tion of degrees of freedom is not always so obvious, since in an itinerant picture all electronic 
bands contribute to both conduction and magnetism with different weights. Thus, the for- 
mulation of spin torques based on the s-d picture is an approximation, and this is a serious 
problem when one tries to evaluate current-induced torques in actual materials. For trust- 
ful estimates, formulations beyond the simple s-d separation is certainly required. Such a 
formalism can be combined with first-principles calculations without any artificial assump- 
tion and would be useful for realistic estimates of current-induced torques and efficiency of 
current-induced switching. The aim of this paper is to develop a new calculational scheme 
satisfying these requirements based on the spin continuity equation. 

Theoretical determination of current- induced torques is difficult even in the simplest case 
of s-d model when spin relaxation and non-adiabaticity is present^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 1 ^. So far, 
very few studies on the effect of spin relaxation due to spin flip scattering by magnetic 
impurities have been done microscopically^ 1 ^. In the s-d formalism, the current- induced 
torque is represented as the effective field due to the spin polarization of the conduction 
electron, s. The torque is therefore given as r^ sd ' = —J sc iS x s, where S is the localized 
(d) electron spin and J s d is the exchange interaction constant. Microscopic calculation using 
linear response theory^ revealed, in agreement with phenomenolocigal result^, that spin-flip 
interaction of conduction electrons with random impurities induces a torque perpendicular 
to the spin-transfer torque (called (5 terms^) . The torque is written as 



where P is the spin polarization of the current and j is the current density. The coefficient 
j3 was calculated by summing over not a few Feynman diagrams representing self-energy and 
vertex corrections^^. 

The case of itinerant ferromagnetism was studied by Tserkovnyak et al— and Duine 
et alM. They introduced the magnetization as a mean-field expectation value of itinerant 
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electron spin, and thus the models considered were effectivelly the s-d model. Tserkovnyak 
et al. considered a kinetic equation for the spin density with a consistency condition for the 
magnetization, but the spin dephasing term was introduced phenomenologically. Duine et al. 
estimated the torques by calculating the effective action for the magnetization fluctuation, 
which has been assumed to be of small amplitude. Within the mean-field treatment, the 
toruqe in the itinerant case turned out to be exactly the same as that of the s-d model^^. 

It has recently been noticed that the coefficient (3 is very important for the realization of 
highly efficient magnetization switching by the current^^^ 1 ^. First, it affects the threshold 
current, and the intrinsic pinning threshold is replaced by an extrinsic one, which is usually 
lower than the intrinsic one. Second, it results in a terminal speed of the wall, v oc "j, 
which can exceed the pure spin-transfer speed limit if £ is large (a is Gilbert damping 
parameter). Third, the deformation of the wall depends on (3. When f3 ~ a, deformation is 
suppressed and weak dissipation may be expected^. Experimental studies of the value of j3 
have recently been carried out. Significant wall deformation observed in permalloy indicated 
that (3 o^. Thomas et a/.— found for permalloy that the observed wall speed corresponds 
to (3 ~ 8a. Therefore, determination of f3 is of particular importance for device applications. 

In this paper, we will present a microscopic calculation scheme different from the s-d 
formalism^. The idea is simply to use the continuity equation of spin, and thus the formu- 
lation is not necessarily based on the s-d interaction picture. The formalism turns out to be 
quite powerful in particular, for the determination of spin relaxation effect, (3. The continu- 
ity equation which we consider is essentially the kinetic equation discussed by Tserkovnyak 
et a/.—, but all observables have been microscopically defined and can be calculated using 
our formalism. For instance, spin dephasing time introduced phenomenolocigally in Ref.— 
is represented by the spin source term (T) defined by Green's function in our formalism. 
Microscopic details of this term T turn out to be essential in determining the spin-relaxation 
induced torque. 

Our scheme is applicable also to the s-d model or mean-field approximation of itinerant 
ferromagnetism. We will use our formalism to obtain the analytical expression of the torques 
arising from both spin-flip scattering and spin-orbit interaction arising from the impurities 
in the s-d model in the adiabatic limit. In the present formalism, the number of contributing 
diagrams is less than the number of diagrams used in the s-d exchange formalism^, and 
thus the calculation is easier. 
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II. FORMALISM 



The spin density s of the total system is defined as the expectation value of conduction 
electron spin, summed over all bands n as 

s a (x,t) = J2(4Xx,t)a a c n (x,t)). (2) 

n 

It satisfies the equation, hs a = i(([H, (4]<7 a c n ) + (c^a^H, c n ])), where H is total Hamil- 
tonian. We assume that H consists of free part, spin relaxation part, H sr , as H = 
J d 3 x ^\^ c n\ 2 + H ST . Then the continuity equation is obtained as 

hs a = --V ■ j s Q + T a , (3) 

where e represents the electron charge. Here the spin current j s is defined by the free part 

as 

Jsp = -^J2( C nM ^ ° a Cn(x',t)) , (4) 
n 

and the spin source (or sink) T is a contribution arising from spin relaxation and interaction, 
i.e., 

T a = zJ2(([H^ci]a a c n ) + (cla a [H SI ,c n })). (5) 

n 

The continuity equation Q is sufficient to calculate the torque acting on the spin. Actually, 
the equation is equivalent to the equation of motion of spin, hs = r, where r represents the 
total torque acting on the spin. The torque is thus simply given by 

r a = --V- j? + T a . (6) 

Note that the continuity equation describes the time-dependence of the spin density, and 
therefore the right-hand side of Eqs. ([3]) and ([6]) is uniquely defined even in the presence 
of spin relaxation, where the spin current can be defined in several different ways, see 
Ref.— . In the context of spin Hall effect, the continuity equation (j3]) was used to obtain 
proper definition of spin current and to explore transport properties^ 1 ^ 1 ^. Concerning the 
current-induced torques, the equation (jHJ) has been so far applied only in the absence of spin 
relaxation term, where the torque is given by the divergence of the spin current^ 1 ^. The 
main aim of this paper is to study the spin relaxation contribution, T. 
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Let us look explicitly at the continuity equation in case of spin relaxation due to spin 
impurities and spin-orbit interaction, H sr = H s f + H so . Spin flip interaction is described by 

H si = v s / d 3 x^ j S imp (x) ■ (4<rc n ), (7) 

where v s is a constant, Si mp (x) = X]™ imp Si mpi S(x — R4), S{ mpi represents the impurity spin 
at x — Ri and n imp denotes the number of impurity spins. The spin-orbit interaction is 
written as 

H so = ~^A so J d 3 x ]T ^jk^jV^\x)(4 l a l Vfe c n ), (8) 

ijkl 

where the potential Vs^ is here assumed to arise from random impurities and depends on 
the spin direction (I). 

The spin-relaxation torque is given by a sum of contributions from spin-flip and spin-orbit 
interaction as T a = T s f + T s ", where 

T s f(x) = 2v s J2e aP7 (sL p si) i (9) 

Pi 

T s a D (x) = -2m\ so MtW W 7) (*)4 • (10) 

The average over random impurity spins and spin-orbit potential is represented by ( ) ; . 

All the terms on the right-hand side of torque Eqs. ([6]), (Q and f flOl) are written in terms 
of local spin density and local spin current, and so the torque acting on the spin is calculated 
by estimating the spin density and the spin current. This representation of the spin torque 
applies to any spin relaxation processes and interaction, and is directly calculable without 
assuming separation of spin and charge degrees of freedoms. Equations (jSJ), © and ( jTUl) are 
thus a suitable starting point for realistic estimates based on first-principles calculations. 
This is the essential point of this paper (although, ab initio calculations using the present 
formalism still need to be undertaken). 

III. APPLICATION TO THE s-d MODEL 

In the latter part of the paper, we will apply this formulation to estimate the current- 
induced torques in the adiabatic limit (i.e., slowly varying magnetization compared with 
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conduction electron motion) to show the validity and usefulness of our formalism. We will 
calculate the torque arising from the spin relaxation due to both the spin flip scattering and 
the spin-orbit interaction. It is found that the torque is represented by the so called the /3 
term in both cases and values of corresponding j3 are calculated. Out formulation is thus 
useful for both the analytical and the numerical studies. 

We will now consider s-d model with only one conduction band. (Please note that the 
assumption of separation of s and d electron here is simply for analytical demonstration and 
is not a requirement of present formulation. ) The s-d interaction between a localized spin 
S and the conduction electrons is given by 



We describe the adiabatic limit by the standard local gauge transformation in the spin 
space, choosing the electron spin quantization axis along S(x,t) at each point. A new 
electron operator a = (a + ,a_) t (t denotes transpose) is defined as c(x,t) = U(x,t)a(x,t), 
where U is a 2 x 2 matrix which we further define as U(x,t) = m ■ cr, m being a real 



as = (m x d il m) a . Then the Hamiltonian of a-electrons is given by the free part, 
J2ka €k c a ka a k(T {^ka = £fc — crM, <j = ± represents the spin), Ha, describing the interaction 
with the SU(2) gauge field, and H em , the interaction with the external electric field which 
drives the currenti^. Here, we consider static local spins in the adiabatic limit, where 
the momentum transfered by the gauge field to conduction electrons is negligibly small 
(compared to hp), and take into account the gauge field only in linear order. Then, the 
gauge interaction is given by^ 




(11) 



three-component unit vector m = (sin | 



| cos (f), sin | sin </>, cos ~ 



|). The gauge field is written 




(12) 



The applied electric field is represented by the interaction 




iehE, 



(13) 



k L ctqr 



where f2o is the frequency of the field chosen as Sl —> at the end of calculation. 
The spin current part of torque is calculated in the adiabatic limit as 



(V^n)j s/i . 



(14) 
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Here, n = S/S and, therefore, this contribution corresponds to the standard spin transfer 
torque. 



A. Torque from spin-flip scattering 

Let us turn to the spin relaxation part of the torque arising from spin impurities, i.e., 
Eq. 0. (The effect of spin relaxation on the spin-current part can be shown to be simply 
due to modification of lifetime, r.) Here, we assume that the impurity spins are influenced 
by a strong s-d exchange field and write S- irap a {x) = R a /3(x)S^ mp (x), where represents 
impurity spin in the rotated frame, and 

R a p = 2m a mp - 5 a/3 , (15) 

is a rotation matrix. Then the averaging is given by [S? (x)^^ p {x')) = \5 a p5{x — 



a O n im P S'imp , where n imp is the impurity spin concentration. (Averaging taken with respect 
to Si m p turns out to lead to essentially the same result as in the case of Si mp .) The spin 
source term is written as 

T s f(x) = —2iv s ^ F a/3j (x) (S£„(x)ta[^G<j) (16) 



Pi 

where 

F a /3^ = ^ ^ £ct[iv Ry^f j (^-^) 

and G^x' = i {o) (x')a(x)^} is the lesser component of the Green's function defined on 
Keldysh contour in the complex time. To the lowest (second) order in v B , we obtain af- 
ter averaging over spin impurities 

T s f(x) = -zjj Ws 2 ^^^F^^ (18) 

where denotes Green's functions without impurity spins but including the gauge field 
A and external electric field E. Including these fields in linear order, we obtain 

T s f(x) = -^n imp v s 2 S~? Yl F aPl {x)E,Al(x)DH\ (19) 

oTYl 

Pltw 



7 



where 



1 f did 

D 1l & = n lim n ?r / tr K { S ^9k'^(r J g kul a s g kuJ+ n 

n -»o "0 7 ^ — 



fcfc' 



i A 4 f 7 5 7 <5 1 A 

H \9k'ujCr gku)9kui+n O' 9ku+Q, + gk'u)& 9kw& gku)9kuj+n ( 

m J 



+ C.C.. 



(20) 



Here, the Green's function g kU ) is the Fourier representation of free Green's function and 
[ ] < denotes the lesser component. They are diagonal in spin space, being defined in gauge- 
transformed space. (Complex conjugates are denoted by c.c.) Figure Q] shows the contri- 
butions to D 1 ^, 6 diagrammatically. The lesser component is calculated in standard manner 



k'uj °L 




k,uj + : jp 




k, uj + ^Vs^4<5 



FIG. 1: The diagrammatic representation of Dffi. Double dashed, dotted and wavy lines denote 
interaction with impurity spin, applied electric field E and gauge field A, respectively. 



in the limit of Qq — * 0. The first two diagrams of Fig. [T] are simplified by use of partial 
integration over k using ^(g k ) 2 = ~ij~9k e ^ c - These contributions are obtained as 



D"f l - 2) = hm /-Vti 



fcfc' 



k,,k 



m 



+ 5 



/lis 

1 

n~ 



{(^(AJV - ^VJV)^ - c.c.} 



/ ( « - Y ) (^fc'^ 7 + 



^0 



(21) 



8 



where fiyj) = (e^ w + 1) 1 . Similarly, the third contribution in Fig. [I] is obtained as 
vM3) - C duJ 



D 



fcfc' 



+ 



{(^G?fc'J V - ^(^J 2 ^)^LAL - c.c. 



+■ ^ (/ - y) (<A^ 7 + ^fc^kL^L 



(22) 



Noting that only antisymmetric part with respect to (3 and 7 contribute to the torque, these 
contributions are summed to be 



DM = _ % 
<* v 2n 



/'Htr [(^Im(^,)^ 7 " ^Mgl'W 



x 



h h 



(kLI 2 A£ + ^l^l 2 ) + Wri Afc 



(23) 



where g k = g r kuJ=0 etc. We see that spin flip processes contributes as additional lifetime as 
indicated by the imaginary part of spin scattered electron Green's function, Img k ,. 

To estimete the trace in the spin space, we use general identities which hold for 2x2 
diagonal matrices B, C, and D (containing only a z and the identity matrix): 



tr^Ba^ - a~<Ba p ){Ca & D + Da s C)] 



tr[(a /3 J Ba 7 - a'Ba^iC^D - Da 6 C)} 



x((BC) + D_ + {BC)-D + + (BD)+C- 
+2e^ z 5 Sz (B + (CD)_ + B.(CD)+)] 
2(<5 72 5/35 — 5p z 5 y s) 

x{(BC) + D^ - {BC)^D + - (BD) + C^ 



{BD)^C + ) 



(BD)_C + ), 
(24) 



where the components B± are defined as B = (B + + £>_ + (B + — B^)a z )/2, etc. The result 



for D^f is then obtained as 



— <W (a(e/3 7 5 - e/3 lz 5 Sz ) + b(5/3s5 JZ - 5 7 s5p z )) 



(25) 
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where the coefficients are given by 



a = — ■ 



-YY 

kk' aa' 



3m 



9kcr\ \9k,-a + 9k -a) + (9k a9k-a + 9ka9k-a) 



( l ™9k>*'] 



b = ~^Yl J2^9ha9k,-a( l ^9k'a'] 



(26) 



kk' aa' 



Using F a/3l = -e al3l -2 J2s m s {e ai sm(3 - e a/ jjm 7 ), and A M = ~(n x d^ri) - A^n^-, the torque 
due to spin-flip is obtained as 



2c 



%i = -— v s 2 S imp ~ ^^(a(n x ^n) - bd^ri). 



(27) 



The coefficients a and 6 are calculated as a = ir(m/e 2 M)(a + — cr_)(z/ + + vJ) and 6 = 
0(a x (e^r)" 1 ) ~ 0, where v± and cr^ = (P'ri^T^jm are the spin-resolved conductivity 
and density of states. (Coefficient b is treated as zero within the present approximation.) 
Therefore, the torque induced by the spin relaxation is simply a (3 term given by 



P 

%t = -Arf — (n x (j ■ V)n), 

e 

where P = (cr + — cr_)/(cr + + a J) is the spin polarization of the current and 

2tt 



(28) 



sf 



3M 



(29) 
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Defining the spin-flip lifetime (r s of Ref.-) as (note that S z + S± of Ref.- corresponds to 

2 2 

fSimp here) r sf _1 = (in / 3)n imp v s 2 S imp (u + + we find /3 sf = ft/( 2iWr sf), which agrees 
with results obtained in Refs.-^. 



B. Torque from spin-orbit interaction 

The torque from spin-orbit interaction, Eq. ( jTUl) . is calculated in a similar way. The 
spin-orbit interaction is written in the rotated frame as 

H so = A so J d 3 x £ijkVjV®(x)R u (x) \^~aW \/ k a + A l k cJc^j . (30) 
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The spin-orbit contributions to the spin current and the electron density in the rotated frame 
are obtained as 



~ p 



x) = -^-Aso^e^V* - V x ')„ j d 3 x 1 V ] V s ^(x 1 )R u (x 1 



x tr 



ijkl 
y(0) 



[ - o(V - V V + A{( Xl ) ) G^ x , 



Ko^tijk, I (fx^jV^^Ra^)^ 1 - V%tr 

ijkl 



+ 0(A). 



(31) 



The torque is then calculated as 



T£{x) = -i\ 2 so ^2 X/ W^j* / d^ii^a^M J] J] m . e -ip-(«-« B1 ) e -i(fc-fcO-» e -i(fci-«4)- a »i 



filiVT jklm 
X (F S o (/3) (P)Ko (T) (-P 



fcfc'p fcifc^ 



(* + k') m Al( Xl )tT Uc^G^J < + (An + A4)*^(*i)tr [^V^,* 



(32) 



In the adiabatic limit we consider, Green's functions are diagonal in wave vectors, G^ k , = 
8k,k'G k \ and the integration over x\ can be carried out treating the slowly varying variables 
R{x{) and A(x\) as constants, resulting in J dx 1 e~ l ( p ~ k+k '* ) '( x ~ Xl ' ) = V5 Pi k~k'- We therefore 
obtain 



T£(x) = -i\ 2 so X X e apLT ei mT e vjk R^{x)R ul {x) ^(k - k')i(k - k') 

fipvr jklm 



kk' 



-{k + k') m {ki + k[)k^ 



x (vJ"\k-k')Vj T \-k + k f 

+ {k + k') m Al{x)tr \<r p G^G^ " + ( k + k') k Ai{x)tr 



(33) 



We average over spin-orbit impurities so that average remains finite only when the spin 
polarizations are parallel. Impurity averaging is thus given as 



Vj v \p)Vj T \-p'))_ = n s J UT S p ^ 



(34) 
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The result of the torque is 



X 



£ 

fcfc' 



(Jfe x k') T (k x fc%tr 



-{k x Jfe') r [(Jfe - Jfe') x A 7 ]^tr 



^(0)^(0)- 

^G^G^] ^ + (k x - k') x A\tr 



(35) 



The last two terms lead to vanishing contribution in the adiabatic limit. In fact, these are 
already linear in A and so does not contain spin-flip conponents, and thus a z and G^ 
commute each other. We therefore obtain 



[(Jfe - k') x A 7 ]„tr 



+ [(Jfe - Jfe') x A%tr 



(5^ z [(k - k') x A\ + 5 ltZ [(k - k') x A^)tr 



a'G^Gf 



(36) 



This contribution is symmetric with respect to (3 and 7, and results in zero when multiplied 
by F^ 1 , which is asymmetric with respect to (3 and 7. 

The first term of Eq. ( 1331) can be simplified by using the rotational symmetry of elec- 
tron, ((Jfe x k') T {k x k') T ) = § ((Jfe x Jfe') • (Jfe x Jfe')) = i ((k 2 k f2 - (Jfe ■ Jfe') 2 )) (( ) denotes the 
angular average), as 



/3/j.ut fcfc' 



(37) 



We therefore see that the expression is similar to that of spin-flip impurity case, Eq. ( fl8l) . 
Including the effect of electric field and gauge field to linear oder in both similarly to the 
spin flip impurity case, we obtain the torque as 



%o = -^n S0 X 80 2 a' E^n x d^n) 



(38) 



where coefficient is given as 



fcfc' 



3m 



9ka\ idk-a + 9k-a) + (dkcrdk-a + 9ko9k-o) 



(39) 
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The coefficient is calculated as a' = tt 3 ^™ M (a + k F+ 2 — <r_k F _ 2 ){y + k F+ 2 + u_k F _ 2 ). There- 
fore, spin-orbit interaction yields the j3 term with coefficient given by 

B =- - (l5±Z± _ lizZh i 

Pso 2Mn + r + -n_r_ I r (») r 0») 1 ' 



where 



27T 

n so \ 2 so k 2 F± (v + k F+ 2 + u_k F J), (41) 



r M 9 

with t±°^ as the lifetime due to spin-orbit interaction. 

The total current-induced torque in the adiabatic limit is therefore given by Eqs. (j!4p 

d2SD (J2HD (SOD as 

r = -£(V ■ j)n - [3 SI j{n x (j • V», (42) 

with (3 SI = (3 si + Pbq. 



IV. SUMMARY 



In summary, we demonstrated that the spin continuity equation represents the current- 
induced torques acting on the magnetization, and that it can be used for microscopic deter- 
mination of the torques. The present formalism does not assume separation of magnetization 
and conduction electron degrees of freedom and can directly be applied to itinerant electron 
systems without mean-field approximation. In this paper, the formalism was applied to the 
s-d model in the presence of spin relaxation caused due to spin-flip scattering and spin-orbit 
interaction with impurities. Both relaxation processes were shown to induce the so called j3 
torque term. 

Application of the formalism to realistic itinerant system using first principles calculations 
would be very interesting, since it would allow for quantitative estimations of current- induced 
switching. Of particular interest are the systems with enhanced spin-orbit interaction near 
surfaces and multilayers. Our formulation can be easliy extended to describe these systems. 

Further improvement of the present theory would be to include effects caused by electron- 
electron correlation. If the correlation is represented within the mean-field approximation 
by a local spin-dependent potential, the torque is straightforwardly calculated similarly to 
the estimate of spin-flip scattering. Treatment beyond mean-field would be an important 
furture work. 
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